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. Abstract 

, Wc report the detection of bursts of 6.7 GHz methanol maser emission in a high-mass star-forming 

region, G33.64— 0.21. One of the spectral components of the maser in this source changed its flux density 
by 7 times that of the previous day, and it decayed with a timescale of 5 days. The burst occurred repeatedly 
I ' in the spectral component, and no other components showed such variability. A VLBI observation with 

2 \ the Japanese VLBI Network (JVN) showed that the burst location was at the southwest edge of a spatial 

distribution, and the bursting phenomenon occurred in a region much smaller than 70 AU. We suggest an 
' impulsive energy release like a stellar flare as a possible mechanism for the burst. These results imply that 

6.7 GHz methanol masers could be a useful new probe for studying bursting activity in the process of star 
formation of high-mass YSOs with a high-resolution of AU scale. 
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. high-mass star-forming regions (Minier et al. 2003; Xu 

• ' A rapid rise of flux called a flare or a burst is often ob- et al. 2008). These masers are thought to be in the 

' served in young stellar objects (YSOs). An X-ray flare regions of a gas density of 10^~® cm~'^ and a dust (in- 
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Introduction methanol masers at 6.7 GHz are observed only around 



with high temperature gas (e.g., Skinner et al. 1997; frared excitation source) temperature of >100 K, corre- 

Tsuboi et al. 1998) and a radio flare with cyclotron or sponding to a distance of <1000 AU from a central star 

gyro-synchrotron radiation (e.g., Phillips et al. 1996; (Cragg et al. 2005). Because the size of a maser spot is 

Massi et al. 2006) have been observed in V773 Tau, as compact as several AU (Minier et al. 2002; Sugiyama 



which is a binary T Tauri (pre-main sequence) star. The et al. 2008a), the maser can be observed using a Very 
\^ ' timescale of these flares is short, on the order of approx- Long Baseline Interferometer (VLBI). The spatial distri- 
5^ , imately 0.1 to 10 days. Such flares are thought to be bution and the radial velocity structure of the maser spots 
caused by magnetic reconnection and its energy release observed by VLBI often show that a disk or a rotating 
around the star (Bower et al. 2003). Because the in- toroid often exists around high-mass YSOs (Minier et al. 
teraction of the gas and the magnetic field at the stellar 2000; Bartkiewicz et al. 2009; Moscadelli et al. 2007), 
surface and the circumstellar environment (disk) is im- and the large part of the methanol maser emanates from 
portant in the process of star formation, it is necessary the gas rotating around or infalling to these YSOs. The 
to reveal the nature of these flares, i.e. magnetic energy flux variability of this maser has been observed with long- 
release, with AU scale spatial resolution. Flare activi- term monitoring observations, and some types of varia- 
ties have been reported mainly for low-mass YSOs so far. tion, such as periodicity, have been found (Goedhart et 
Although some authors (Tsuboi et al. 1999; Hamaguchi al. 2003, 2004, 2009). In some sources observed by flux 
et al. 2005; Stelzer et al. 2005) discussed that some part monitoring and spatial distribution with VLBI, it is pro- 
of X-ray variability were magnetic origin, high- mass YSOs posed that the variability of the maser is caused by the 
are not considered active in terms of bursting variability. variation of the exciting source (Sugiyama et al. 2008b; 

Wc studied the 6.7 GHz methanol maser which helped Araya et al. 2010; Szymczak et al. 2011). 
us understand not only the gas dynamics but also the When a stellar flare occurs, the released energy would 

possible flaring activity around high-mass YSOs. The heat the gas and dust around the flare region, cans- 
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ing a maser flare. The maser flare can be used to re- 
veal the location at which magnetic reconnection and/or 
other activity occurs around the YSO by conducting a 
VLBI observation. To look for such maser flares, we se- 
lected 22 maser sources that have relatively large flux 
densities and multiple spectral components, and carried 
out a daily monitoring observation. The detail of the 
whole observation will be reported elsewhere, we report 
here focusing on the discovery of a bursting methanol 
maser in G33.64-0.21 (IRAS 18509-^0027, 18'^53™32f563 
-h00°31'39"18, J2000.0 after Bartkiewicz et al. 2009). 

This source is a high-mass star-forming region located 
at the near-kinematic distance of 4.0 kpc, and the in- 
frared luminosity obtained from the IRAS database was 
1.2 X lO''^ Lq. Although there are some infrared sources 
detected with 2MASS, Sptzer/GLIMPSE and MSX near 
the methanol position, it is not clear which one is the as- 
sociating source due to the large uncertainty of the maser 
position {'^ 5"). The methanol maser of this source was 
first observed by the Torun survey (Szymczak et al. 2000). 
There are at least five spectral features (components I-V, 
defined in radial velocity order) with a typical line width 
of 0.3 km and one weak, broad feature (component 
VI) in the spectrum (figure 1). The magnetic field of this 
source was measured through polarization observation as 
— 18 niG, which is stable over the entire velocity range 
(Vlemmings 2008). 

In this paper, we describe the spectra monitoring and 
VLBI observations and the results in section 2, then we 
discuss interpretations of the mechanism of the bursting 
methanol maser in section 3. 

2. Observations and Results 

2.1. Spectra Monitoring 

The monitoring observations were made with the 
Yamaguchi 32 m radio telescope from July 2 (Day-of- 
Year 183) to October 17 (DOY 290), 2009. A total of 
62 observations were made over the 108-day observation 
period. The observations were usually made every one or 
two days, and the maximum interval was five days. The 
spectrometer used for this observation was the IP- VLBI 
system (Kondo et al. 2003). The sampling was made at a 
rate of 8 Msamples sec~^ (bandwidth 4 MHz) and at two 
bits per sample. The spectra were obtained by a software 
spectrometer with 16384 spectral channels, then binned 
into 4096, yielding a velocity resolution of 0.044 km s~^. 
The integration time was 14 minutes. The left and right 
circularly polarized waves were observed simultaneously 
and were averaged after being transformed into the spec- 
tra. A Icr rms noise level is 1.2 Jy. The accuracy of the 
absolute value of the flux density is approximately 10%, 
while the standard deviation of the relative fluctuation is 
3% estimated from the daily data. The Vlsr calculation 
may contain <0.3 km constant offset. 

Two rapid rises in the flux density (hereafter, burst) 
were detected in G33.64-0.21 at DOY 186 and 211 only in 
component II (T^lsr = 59.6 km s~^), as shown in figure 2. 
At the first burst, the flux density rose to 7 times within 



- DOY 210 
-D0Y2I1 




I II III IV 



56 57 58 59 60 61 62 63 64 65 



Fig. 1. Spectra of G33.64-0.21. The spectra of G33.64-0.21 
observed with the Yamaguchi 32 m radio telescope are shown. 
The dotted and soUd lines show the spectra of day-of-year 
(DOY) 210 and 211, respectively, just before and after the sec- 
ond burst. Only spectral component II (Vlsr = 59.6 km s~^) 
rose by a factor of seven (26 to 191 Jy). 

three days (DOY 183 to 186), and it fell exponentially 
with an e-folding time of 5 days. Because there was three 
days lack of observation just before the burst, the rising 
timescale and the peak flux density contain some ambigu- 
ity. In the second burst, the flux density rose from 26 Jy 
to 191 Jy, i.e., a factor of 7 times within 24 hours from 
DOY 210 to 211, and then it fell slowly The decrement of 
the flux density was exponential, with an e-folding time of 
5 days. This behavior is similar to that of the first burst. 
A small burst was observed at DOY 262. There was no 
change whatsoever in the other spectral components when 
the bursts occurred. The flux density of component I de- 
creased from 75 to 20 Jy from DOY 208 to DOY 240. Such 
variation with a timescale of 20 days has been observed 
in other sources (Sugiyama et al. 2008b; Goedhart et al. 
2009). 

The peak velocity of component II was measured 
by gaussian fitting. The velocity slightly decreased 
(0.05 km s~^) when the component bursted, and the ve- 
locity returned to its original value after the burst (figure 
3). This velocity change repeated for the two bursts with 
an accuracy of less than 0.01 km s~^. This velocity shift 
may be explained as blending two spectral components 
with almost the same but slightly different velocities; one 
component has high velocity and stable flux, and the other 
has low velocity and shows the bursting variability. This 
small velocity variation and its repeatability suggest that 
the bursting region is not disturbed by violent gas motion 
like a Shockwave even when they showed large bursts. 

Szymczak et al. (2010) reported that G33.64-0.21 
showed a large flare only in the spectral component II 
(59.6 km s~i) near JD=2455160 (DOY 330 of 2009). The 
observed ratio of the maximum and minimum of the flux 
density was 16.7. This suggests that G33.64— 0.21 bursted 
again at this time. They also reported the decay time of 
'~ 30 days, however, this might be an wrong estimation 
due to their sparse observation (intervals of 1 to 6 weeks). 
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Fig. 2. The variation of spectral eomponents. The changes 
of peak flux density for each six spectral components (I- VI) 
are shown. Only component II (filled circle and the solid line) 
shows the burst, while the other components are completely 
unrelated to the burst. Two large bursts (DOY 186, 211) and 
a small burst on DOY 262 were seen during the 108 days of 
observation run. 



Fig. 3. Peak velocity and flux density of component II. The 
change of peak velocity of component II with respect to its 
flux density is shown. A negative correlation is obvious. Filled 
circles are for the first burst (DOY 183 to 210), and open cir- 
cles arc for the second burst (DOY 211 to 248), respectively. 
The velocity changes repeated for the two bursts with an ac- 
curacy of less than 0.01 km s~^. 



A bursting phenomenon hke the one observed in 
G33.64— 0.21 was found only in this source among the 
22 sources observed. A large bursting variability in 
methanol masers was reported for NGC6334F (Goedhart 
et al. 2004). The variation pattern of NGC 6334F, i.e., 
the fast rise and relatively slow fall, looks like the burst 
of G33.64— 0.21; however, the timescale of the burst in 
NGC6334F is as long as 100 days. Periodic flare phenom- 
ena have been found in some sources of methanol maser. 
Araya et al. (2010) reported a periodic flare of formalde- 
hyde and methanol masers. The duration of the flare is 
- 30 days for IRAS 18566-K0408, while it is only 5 days for 
G33.64— 0.21. Szymczak et al. (2011) discussed the time 
profile of the flares of two sources; G22. 357-1-0. 066 and 
G9.62+0.20E (the data is from Goedhart et al. 2004). 
The rise and fall (decay) time are 14.5 ± 1.7 and 42.9 ±0.7 
days for G22. 357-1-0. 066. These are ten times larger than 
those of G33.64— 0.21. The ratios of the rise and fall time 
are 0.34 and 0.37 for G22.357-I-0.066 and G9.62+0.20E, re- 
spectively, while it is less than 0.2 for G33.64— 0.21. The 
time asymmetry is much obvious for G33.64— 0.21. 

2.2. VLB I Observation 

A VLBI observation using the Japanese VLBI Network 
(JVN) was conducted on October 6, 2009 (DOY 279) 
from 05:00 to 14:00 UT, to determine the position in 
G33.64— 0.21 where the burst occurred. Although the 
burst did not occur at this time, the bursted component 
(II) was visible. Five JVN telescopes (Yamaguchi 32 m 
and four VERA 20 m telescopes) were used. The obser- 
vation frequency was 6665-6669 MHz, and the number of 
frequency channels was 1024. The spatial distribution of 
the 35 maser spots is shown in figure 4. Those maser spots 
were detected with a signal-to-noise ratio larger than 5, 
and at least 2 adjacent frequency channels. The synthe- 
sized beam size was 11.7x3.0 mas^ with a position angle 
of —45.2 deg, as shown in the lower-left corner of the fig- 
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Fig. 4. Distribution of maser spots observed with the JVN. 
The spot size and color indicate its peak intensity in loga- 
rithmic scale and its radial velocity (see color index at the 
right), respectively. The origin of this figure is the position 
of the strongest component at Vlsr = 60.23 km s~^. The 
maser spots are distributed in an arched structure with an 
extension of 650x200 AU. Labels of I to VI correspond to the 
spectral components. The maser spots of the bursting spec- 
tral component (II) are located at the southwest edge of the 
distribution, as indicated with the dashed ellipse. 

ure. The rms noise of the image in a line-free channel is 
70 mJy beam~^. 

The maser spots form an arched structure with a size 
of 160 mas, corresponding to the linear size of 650 AU at 
4.0 kpc. This distribution is essentially the same as the 
result by Bartkiewicz et al. (2009). There is no peculiar 
characteristic in the spatial distribution of G33.64— 0.21 
compared with those of other sources shown by VLBI sur- 
veys (Minier et al. 2000; Dodson et al. 2004; Sugiyama 
et al. 2008a; Bartkiewicz et al. 2009). Component II, 
which showed the burst, is located at the southwest edge 
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of the spatial distribution of the maser spots. There are 
two clusters of spots that are 30 AU apart in this region 
comprising spectral component II: it is unknown which of 
them caused the burst. The VLBI image reconstructed 
the flux density of 62% of the Yamaguchi spectrum for 
the component 11. This suggest that the size of the com- 
ponent II is as small as the beam size (11.7x3.0 mas^). 
The spots of spectral components I and IV are located 
70 AU in northeast and 130 AU in east direction away 
from the clusters of component II. Because spectral com- 
ponents I and IV did not show any change when compo- 
nent II bursted, the burst-related phenomenon occurred 
within a spatial scale of much less than 70 AU. 

3. Discussions 

The observed properties of the burst are summarized 
as follows: (1) the rising timescale was less than one day, 
whereas the falling was 5 days; (2) the burst only oc- 
curred in spectral component II, at the southwest edge of 
the spatial distribution; and (3) the bursting phenomenon 
occurred in a region much smaller than 70 AU. Before con- 
cluding that the magnetic energy release was the plausible 
origin of the burst, other possibilities are considered. 

It was reported that interstellar scintillation caused the 
variability of a hydroxyl mascr (Clegg & Cordes 1991). 
The amplitude (7 times) and time scale (less than 1 day) 
of the burst in G33.64— 0.21 are, however, the strong ar- 
guments against that the scintillation is the origin of the 
burst. In addition, if the scintillation was the origin of the 
burst, all spectral components would show some variabil- 
ity; the bursts were observed only in one spectral compo- 
nent of G33.64— 0.21. Therefore, the burst was not caused 
by interstellar scintillation. 

In some maser sources, the flux variability of the 
methanol maser is proposed to be caused by the variability 
of the luminosity of the exciting source (Sugiyama et al. 
2008b; Araya et al. 2010; Szymczak et al. 2011). The vari- 
ation pattern of the flux density should be synchronized or 
correlated between two or more spectral components if the 
maser variation is caused by the luminosity variation of 
the exciting source. In the case of G33.64— 0.21, however, 
not even the spectral component nor maser spots spatially 
adjacent to the bursting spots within a distance of only 
70 AU varied when component II bursted. Therefore, the 
burst of G33.64— 0.21 was not caused by a variation of the 
exciting source. 

Shimoikura et al. (2005) showed by VLBI monitor- 
ing that the burst of the 22 GHz water maser in Orion 
KL in 1998 was caused by an overlapping and amplify- 
ing effect that occurred between two maser clouds, which 
were on the line of sight by chance. If we apply the 
overlapping model to the burst of G33.64— 0.21, its short 
timescale would be a difficulty. Because the spectrum of 
G33.64— 0.21 has a velocity range of 6 km s~^, one can 
assume 3 km s~^, half of the range, to be the speed of 
the gas cloud. The maser cloud can move only 3 x 10^ m 
in one day at such speed, corresponding to 0.5 /las at a 
distance of 4.0 kpc. The brightness temperature would ex- 



ceed 10^^ K if a flux rise of 165 Jy occurred in this small 
region. This brightness temperature cannot be achieved 
with the radiation mechanism of the 6.7 GHz methanol 
maser. Moreover, the asymmetry of the variation of rise 
(timescales of less than 1 day) and fall (5 days) cannot be 
explained with this model. 

The simplest explanation of this burst is that energy 
was injected within a short time into the small region at 
which the burst occurred. An example of such energy re- 
lease seen in astrophysics is the solar flare, which is an 
impulsive event of magnetic fleld reconnection. The en- 
ergy accumulated in the magnetic field on the solar sur- 
face is discharging by the reconnection of the magnetic 
field in a short time. The released energy causes a parti- 
cle acceleration that heats the gas near the reconnection 
position, and finally, radiation is emanated from the gas. 
Similar to the solar flare, flares at a much larger scale 
are observed in T Tauri stars, such as V773 Tau (Phillips 
et al. 1996) or GMR-A (Bower et al. 2003; Furuya et 
al. 2003). Here, we propose a mechanism of the maser 
burst: an impulsive energy release which resembles those 
observed in active YSOs occurred at the stellar surface 
or in the circumstellar disk around the high-mass YSO 
in G33.64— 0.21, causing the maser burst. The released 
energy was converted into the heat of the gas and dust 
surrounding the release point. The population inversion 
of the methanol molecule was strengthened by the infrared 
emission from the heated dust. Finally, the maser emis- 
sion emanating from the gas around the release point was 
suddenly strengthened. The asymmetry of the rise and 
fall of the burst is naturally explained within this model 
because the rise is impulsive, but the fall depends on the 
relatively slow radiative cooling process. The cooling time 
scale was theoretically estimated of 1.2 days by van der 
Walt et al. (2009) for the optically thick condition. This 
time scale was found to be too short to explain the decay 
time of NGC6334F maser flare 100 days), but could be 
applicable to the case of G33.64— 0.21, which has signifi- 
cantly shorter decay time 5 day) . Because the velocity 
change during the burst was small, the shock process is 
not suitable to explain this burst. A theoretical calcu- 
lation (Cragg et al. 2005) would support our scenario 
because the brightness of the maser depends on the dust 
temperature, while the gas density and gas temperature 
are insensitive parameters to the brightness temperature 
variation. A change of dust temperature of 100 K (e.g. 
150 to 250 K) would cause the increment of maser flux of 
7 times. 

A brief consideration of energy dynamics is presented 
below. Assuming the size of the area where the gas is 
heated is 10 AU in diameter and the number density of 
the hydrogen molecule is 10^ cm""^, then the gas mass in 
this area would be 6 x 10^^ kg. An energy of 6 x 10^*" J 
is required to heat the gas and raise its temperature by 
100 K so that the mascr excitation condition changes dras- 
tically. This energy is comparable to the energy of an 
X-ray flare observed in V773 Tau (10^^ J; Tsuboi et al. 
1998). If we assume that the energy was released from the 
magnetic field in/around the star, and that the magnetic 
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field strength was 18 niG throughout this region, the to- 
tal magnetic energy in this volume would be 2 x 10'^° J. 
Because the magnetic energy is much higher than the re- 
quired energy for the burst, the magnetic energy release 
can account for the burst. 

The large bursts occurred twice, and they only oc- 
curred in spectral component II, which means that the 
energy releases repeatedly occurred at the southwest edge 
of G33.64— 0.21. At the bursting place, there might be 
a binary companion or planet, which could cause a local 
disturbance and trigger the energy release. The distribu- 
tion of a water maser of G33.64— 0.21 (Bartkicwicz et al. 
2011) is elongated north and south from the position near 
the bursting methanol component (II), but the positional 
relation between the water maser and the central star is 
not definitive because of the large positional uncertainties 
of 5" and 0.15" of the methanol and water maser, respec- 
tively. 

If the maser burst was caused by local heating, a 
local burst would be observed simultaneously at in- 
frared and submillimetre bands. Such observation re- 
quires high spatial resolution, which will be achieved by 
ALMA (Atacama Large Millimeter/submillimeter Array). 
Finally, we emphasize that this result demonstrates the 
6.7 GHz methanol maser can be a useful tool to investi- 
gate the bursting activity in the process of star formation 
around high-mass YSOs. 

4. Conclusion 

We have presented the discovery of the bursting activ- 
ity in a high-mass star-forming region detected with the 
6.7 GHz methanol maser. A 108-day continuous moni- 
toring for 22 sources of methanol maser was made with 
the Yamaguchi 32 m radio telescope. In one source, 
G33.64— 0.21, rapid rises of the flux density (burst) were 
observed twice only in one spectral component out of 6 
components. In the second burst, the flux density rose 
7 times within 24 hours, and then it fell exponentially 
with an e-folding time of 5 days. There was absolutely 
no change in the other spectral components when the 
bursts occurred. A VLBI observation using the Japanese 
VLBI Network (JVN) showed that the maser spots form 
an arched structure and the burst component is located 
at the southwest edge of the spatial distribution of the 
maser spots. There are spots which did not show burst 
70 AU from the burst spot. This means that the burst- 
related phenomenon occurred within a spatial scale of 
much less than 70 AU. In order to account these prop- 
erties, we proposed a mechanism of the maser burst as 
an impulsive energy release like the stellar flare observed 
in YSOs occurred at the stellar surface or in the circum- 
stellar disk in G33.64— 0.21. Finally, we noted that the 
6.7 GHz methanol maser can be a useful tool to investi- 
gate the bursting activity in the process of star formation 
around high-mass YSOs. 
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